Astrophysics independent bounds on the annual modulation of dark matter signals 
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We show how constraints on the time integrated event rate from a given dark matter (DM) direct 
detection experiment can be used to set a stringent constraint on the amplitude of the annual 
modulation signal in another experiment. The method requires only very mild assumptions about 
the properties of the local DM distribution: that it is temporally stable on the scale of months and 
spatially homogeneous on the ecliptic. We apply the method to the annual modulation signal in 
DAMA/LIBRA, which we compare to the bounds derived from the constraints on the time-averaged 
rates from XENON10, XENON100, CDMS and SIMPLE. Assuming a DM mass of 10 GeV, we 
show that a DM interpretation of the DAMA/LIBRA signal is excluded at 6.3a (4.6a) for isospin 
conserving (violating) spin-independent interactions, and at 4.9a for spin-dependent interactions on 
protons. 



Introduction. Dark matter (DM) constitutes a sig- 
nificant fraction of the energy density in the universe, 
Hdm = 0.229±0.015 pp. This conclusion is based entirely 
on gravitational effects of DM. A fundamental question is 
whether DM interacts also non-gravitationally. There are 
a number of experiments searching for signs of such DM 
interactions. Direct detection experiments, for instance, 
are looking for a signal of DM particles from the galac- 
tic halo that would scatter in underground detectors. A 
characteristic feature of the resulting signal will be an 
annual modulation, because the earth rotates around the 
sun, while at the same time the sun moves relative to the 
DM halo g]. 

At present two experiments are reporting annually 
modulated signals, DAMA/LIBRA [3] (DAMA for short) 
and CoGeNT [4], with significances of 8.9<r and 2.8a, re- 
spectively. Are these signals due to DM? The answer is 
readily obtained by (i) assuming a specific local DM ve- 
locity distribution and (ii) postulating the predominant 
DM-nucleus interaction. Usually a simple Maxwellian 
DM halo is adopted. If interpreted in terms of elastic 
spin-independent DM scattering both claims are in ten- 
sion |6] with bounds on time integrated rates from 
other direct detection experiments such as XENON10 [7], 
XENON100 |S], or CDMS [5]. The situation may change 
in the case of non-standard DM halos with, e.g., highly 
anisotropic velocity distributions, DM streams or DM de- 
bri flows. Recently CDMS provided a direct bound on the 
modulation signal, which disfavors the CoGeNT modu- 
lation without referring to any halo or particle physics 
model [10 . Therefore we focus below mainly on DAMA. 

In this Letter we present a general method that avoids 
astrophysical uncertainties when comparing putative DM 
modulation signals with the bounds on time averaged 
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DM scattering rates from different experiments. The 
method combines the results from [TT1 [12] with the 
bounds on the modulation derived by us in [13] . We 
are then able to translate the bound on the DM scatter- 
ing rate in one experiment into a bound on the annual 
modulation amplitude in a different experiment. The re- 
sulting bounds present roughly an order of magnitude 
improvement over [LTl [Dl] and |13| . 

The bounds are (almost completely) astrophysics in- 
dependent. Only very mild assumptions about DM halo 
properties are used: (i) that it does not change on the 
time-scales of months, (ii) that the density of DM in the 
halo is constant on the scales of the earth-sun distance, 
and (iii) that the DM velocity distribution is smooth on 
the scale of the earth velocity v e = 29.8 km/s. If the 
modulation signal is due to DM, then the modulation 
amplitude has to obey the bounds. In the derivation an 
expansion in v e over the typical DM velocity ~ 200 km/s 
is used. The validity of the expansion can be checked 
experimentally, by searching for the presence of higher 
harmonics in the time-stamped DM scattering data [12] . 

Bounds on the annual modulation. We focus 
on the case of DM \ elastically scattering off a nucleus 
(A, Z) and depositing the nuclear recoil energy E nr in 
the detector. The differential rate in events/keV/kg/day 
is then given by 

R A (E nr ,t) = Px °\ F\{E nr ) r]{v m ,t) , (1) 

with p x the local DM density, a A the total DM-nucleus 
scattering cross section at zero momentum transfer, m x 
the DM mass, and fi x A the reduced mass of the DM- 
nucleus system. FA{E nr ) is a nuclear form factor. For SI 
interactions with a nucleus (A, Z), a A can be written as 
af = a p [Z +{A- Z)(f n /f p )] 2 ii xA /(i xp , where a p is the 
DM-proton cross-section and /„ iP are coupling strengths 
to neutron and proton, respectively. Apart from a com- 
mon overall factor p x the astrophysics enters the pre- 
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dieted rate in Eq. |T]) through the halo integral 
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where v m is the minimal velocity required to have at 
least E nr energy deposited in the detector. The func- 
tion /(j e t(v,t) describes the distribution of DM particle 
velocities in the detector rest frame with /d e t(v,i) > 
and J d 3 w/det( v 7 1) = 1. It is related to the velocity dis- 
tribution in the rest frame of the sun by /<j e t(v,t) — 
Jsun(v + v e (t)), where v e (t) is the velocity vector of the 
earth. The rotation of the earth around the sun intro- 
duces a time dependence in the DM-nucleus scattering 
rate through r/(v m ,t) = fj(v m ) + Srj(v m ,t), where 



Sr)(v m ,t) = A r ,(v m )cos2ir[t - t (E nr )}, 



(3) 



when expanding to first order in v c = 29.8 km/s <C v sun ~ 
230 km/s. Here, A v (v m ) is defined to be positive. 

Let us now assume that f S \m(v) is smooth on the scale 
of i> e , and the only time dependence comes from the ro- 
tation of the earth around the sun and f S un(v) itself is 
constant in time and space. Then the modulation am- 
plitude A v (v m ) can be bounded in terms of the unmod- 
ulated halo integral fj in the following way |13j : 



A v {v m ) < v e 



drj 



dv 



(4) 



The first term in Q is positive since rj(v m ) is a 
monotonously decreasing function of v m . If we further 
assume that the DM halo is symmetric, so that there 
is only one single direction related to the DM flow (see 
|13| for details), then one obtaines a more stringent con- 
straint: 



dv m A v (v m ) < sin a w e 



rj(vi) - vi dv 



r](v) 
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Here a is the angle between the DM flow and the direc- 
tion orthogonal to the ecliptic. The most conservative 
bound is obtained for sin a — 1 (which would correspond 
to a DM stream parallel to the ecliptic). However, in 
many cases the DM flow will be aligned with the motion 
of the sun within the galaxy. This holds for any isotropic 
velocity distribution and, up to a small correction due to 
the peculiar velocity of the sun, also for tri-axial halos or 
a significant contribution from a possible dark-disc. In 
this case we have sin a ~ 0.5. 

In the following we will use time averaged rates from 
various experiments to derive an upper bound on fj(v m ). 
In order to be able to apply this information we integrate 
Eq. Q over v m and drop the negative terms in Eqs. Q 
and pj). This gives the bounds 
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In practice the integrals on the l.h.s. are replaced by a 
sum over bins. Below we will refer to the relations ([6| and 
Q as the bounds from "general halo" and "symmetric 
halo" (where we will take sin a = 0.5), respectively. 

Bounds on the unmodulated halo integral. Let 
us first consider SI scattering with /„ = f p . Generaliza- 
tion to isospin violating scattering with /„ ^ f p and to 
SD scattering is straightforward. The predicted number 
of events in an interval of observed energies [Ei , E 2 ] is 
given by 

POO 

N \eTe*\= MTA2 / dE nr Fl(E nr )G [El , E2] (E nr )fi(v m ). 
Jo 

(8) 

Here G[E 1 .E 2 ](E n r) is the detector response function, 
which describes the contribution of events with true 
nuclear-recoil energy E nr to the observed energy interval 
[Ei,E 2 ]. It may be non-zero outside the E nr e [E 1 ,E 2 ] 
interval due to the finite energy resolution and includes 
also (possibly energy dependent) efficiencies. M and T 
are the detector mass and exposure time, respectively, 
and we defined 
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where fj has units of events/kg/day/keV. 

Now we can use the fact that fj is a falling function [TT] 
(see also [TU [15] ) . Among all possible forms for fj such 
that they pass through fj(y m ) at v m , the minimal number 
of events is obtained for fj constant and equal to fj(v m ) 
until v m and zero afterwards. Therefore, for a given v m 

](v m ) > n(v m ) with 



we have a lower bound NPl cd ^ 



where E(v m ) is given in (pj). 



fx(v m ) = MTA 2 fj(v m ) / dE nr Fl(E nr )G[ Elj E 2 ](E nr ). 
Jo 

(10) 

Suppose an experiment 
observes N°g* E ^ events in the intervall [Bi,^]- Then 
we can obtain an upper bound on fj for a fixed v m at 
a confidence level CL by requiring that the probability 
of obtaining N°£* E ^ events or less for a Poisson mean 
of /i(u m ) is equal to 1— CL. Note that this is actually a 
lower bound on the CL, since Eq. ( 10 ) provides only a 
lower bound on the true Poisson mean. For the same 
reason we cannot use the commonly applied maximum- 
gap method to derive a bound on fj. If several different 
nuclei are present, there will be a corresponding sum in 
Eqs. pi and (pi). 



The limit on fj can then be used in the r.h.s. of Eq. ^ 
or |7| to constrain the modulation amplitude. For con- 
creteness we first focus on the annual modulation in 
DAMA. If m x is around 10 GeV, then DM particles do 
not have enough energy to produce iodine recoils above 
the DAMA threshold. We can thus assume that the 
DAMA signal is entirely due to the scattering on sodium. 
We define A v = a v p x j (2m x Pp)A ri , which is related to the 
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observed modulation amplitude A° hs by 
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Here <?Na = dE ee /dE nr is the sodium quenching fac- 
tor translating keVee into keVnr, for which we take 
<ZNa = 0.3. The index i labels energy bins, with v l m given 
by the corresponding energy bin center using Eq. @. 
Further, (F^ a )i is the sodium form factor averaged over 
the bin width and /n s = TiNa/^Na + mi) is the sodium 
mass fraction of the Nal crystal. For the modulation am- 
plitude in CoGeNT we proceed analogously. Note that 
the conversion factor from fj to fj is the same as for A n to 
Arj, and does not dependent on the nucleus. Therefore, 
the bounds (|6| and (|7) apply to fj, A v without change, 
even if the l.h.s. and r.h.s. refer to different experiments. 

Let us briefly describe the data we use to derive the up- 
per bounds on fj. We consider results from XENON10 [7 
(XE10) and XENON100 [5] (XE100). In both cases we 
take into account the energy resolution due to Poisson 
fluctuations of single electrons. XE100 is sensitive to the 
interesting region of v m only because of upward fluctu- 
ations from below the threshold. We adopt the best-fit 
light-yield efficiency L c g from [5]. The XEI0 analysis 
is based on the so-called S2 ionization signal which al- 
lows to go to a rather low threshold. We follow [7] and 
impose a sharp cut-off of the efficiency below the thresh- 
old. From CDMS we use results from a dedicated low- 
threshold (LT) analysis [9] of Ge data, as well as data on 
Si I6J. In the case of SD scattering on protons particu- 
larly strong bounds are obtained from experiments with 
a fluorine target. We consider the results from SIMPLE 
[FT] , which uses F5C2CI. We use the observed number 
of events and expected background events to calculate 
the combined Poisson probability for Stage I and 2. For 
the prediction we include energy dependent threshold ef- 
ficiencies from [17]. 

For all experiments we use the lower bound on the 
expected events, Eq. ( [Tpj ) , to calculate the probability of 
obtaining less or equal events than observed. For XE100, 
CDMS Si, and SIMPLE we just use the total number of 
events in the entire reported energy range. For XE10 and 
CDMS LT the limit can be improved if data are binned 
and the corresponding probabilities for each bin are mul- 
tiplied. This assumes that the bins are statistically in- 
dependent, which requires to make bins larger than the 
energy resolution. For XE10 we only use two bins. For 
CDMS LT we combine the 36 bins from Fig. 1 of [9] into 
9 bins of 2 keV where the energy resolution is 0.2 keV. 

Results. In Fig. [I] we show the 3er limits (CL = 
99.73%) on fj compared to the modulation amplitudes 
A\ from DAMA and CoGeNT for a DM mass of 10 GeV. 
Similar results have been presented in [21 [T5] . The Co- 
GeNT amplitude depends on whether the phase is floated 
in the fit or fixed at June 2nd [6], which applies to the 
"general" and "symmetric" halos, respectively. Already 
at this level XE100 is in tension with the modulation 
from DAMA (and to some extent also CoGeNT). 
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FIG. 1: Upper bounds on f\ at 3cr from XENON100, 
XENON10, CDMS LT, CDMS Si, and SIMPLE. The modu- 
lation amplitude A v is shown for DAMA (for <jrjv a = 0.3) and 
CoGeNT for both free phase fit (general) and fixing the phase 
to June 2nd (symmetric). We assume a DM mass of 10 GeV 
and SI interactions. 
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FIG. 2: Integrated modulation signals, f"^ dvA^, from 
DAMA and CoGeNT compared to the 3<r upper bounds for 
the general halo, Eq. We assume SI interactions and a 
DM mass of 10 GeV. The integral runs from v\ = u m i n 
v 2 = 743 km/s (end of the 12th bin in DAMA). 
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We now apply our method. As shown in Fig. [2] the 
null search results become significantly more constrain- 
ing after applying the bounds on the integrated annual 
modulation /J 2 dvA f) from Eq. DAMA and GoGeNT 
are strongly excluded by the bounds from XE100, XE10, 
CDMS LT even for the general halo. If one were to as- 
sume in addition that the halo is symmetric, the bounds 
would get even stronger. Then also CDMS Si excludes 
DAMA, and there is some tension with SIMPLE (not 
shown) . 

In Fig. [3] we consider two variations of DM-nucleus 
interaction. The upper panel is for the case when the 
DM particle couples to the spin of the proton. The null 
search result of Xe and Ge experiments are then irrel- 
evant. However, the bound from SIMPLE is in strong 
disagreement with the modulation signal in DAMA, due 
to the presence of fluorine in their target. (A compa- 
rable limit from fluorine has been published recently 
by PICASSO 18 j.) In the lower panel of Fig. [3] we 
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FIG. 3: Integrated modulation signal P 
DAMA compared to the 3a upper bounds for the general halo, 
Eq. |6f (solid), and symmetric halo, Eq. Q with sin a = 0.5 
(dotted). We assume a DM mass of 10 GeV, and SD in- 
teractions on protons (upper panel) and SI interactions with 
fn/fp = —0.7 (lower panel). The upper limit of the integra- 
tion is V2 = 743 km/s. 



show the case of SI isospin violating interactions with 



fn/ fp 



-0.7. This choice evades bounds from Xe, but 



now the DAMA modulation is excluded by the bounds 
from CDMS Si for the general halo and CDMS Si, LT, 
and SIMPLE for the symmetric halo. 

Let us now quantify the disagreement between the ob- 
served DAMA modulation and the rate from another 
null-result experiment using our bounds. We first fix v m . 
To each value of fj(v m ) Eq. (10 1 provides a Poisson mean 
/j,(v m ). We can then calculate the probability p v to ob- 
tain equal or less events than measured by the null-result 
experiment. Then we construct the bound on the modu- 
lation using the same value fj(v m ) on the r.h.s. of Eq. |6]) 
or (the integrand fj(v) in Eq. ^ is calculated using 
the same p v but with v > v m in Eq. (fl0|)). We calculate 
the probability pa that the bound is not violated by as- 
suming on the l.h.s. of Eq. (|6| or |7| a Gaussian distribu- 
tion for the DAMA modulation signal with the measured 
standard deviations in each bin. Then Pj int(??) = Pr/PA 
is the combined probability of obtaining the experimen- 
tal result for the chosen value of fj. Then we maximize 
Pjoint('7) with respect to fj to obtain the highest possible 
joint probability. 

The results of such an analysis are shown in Fig. [4] 
The analysis is performed at the fixed v m corresponding 
to the 3rd modulation data point in DAMA, depending 
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FIG. 4: The probability that the integrated modulation am- 
plitude in DAMA (summed starting from the 3rd bin) is com- 
patible with the bound derived from the constraints on fj for 
various experiments as a function of the DM mass. The label 
SI (SD), refers to spin-independent (spin-dependent) interac- 
tions with /„ = f p (/„ = 0) , and IV refers to isospin- violating 
SI interactions with fn/fp = —0.7. For solid and dashed 
curves we use the bounds from Eqs. (|6| and (JT]), respectively. 



on the DM mass m x . We find that for all considered in- 
teraction types and m x < 15 GeV at least one experiment 
disfavors a DM interpretation of the DAMA modulation 
at more than 4ct even under the very modest assumptions 
of the "general halo" . In the case of SI interactions the 
tension with XE100 is at more than 6a for m x > 8 GeV 
and saturates at the significance of the modulation data 
point itself at about 6.4<r for m x > 13 GeV. The exclu- 
sion from XE10 is nearly independent of the DM mass 
slightly below 6er. We show also a few examples of the 
joint probability in case of a "symmetric halo" (dashed 
curves) . 

While astrophysics uncertainties are avoided, the ob- 
tained bounds are still subject to nuclear, particle physics 
and experimental uncertainties. For instance, the ten- 
sion between the DAMA signal and the bounds depends 
on the value of the Na quenching factor (jNa, light yield 
or ionization yield efficiencies in Xe, upward fluctuations 
from below threshold, and so on. For example, if a value 
of qNa — 0.45 is adopted instead of the fiducial value of 
0.3 consistency for SD and isospin violating interactions 
can be achieved in the case of the general halo at around 
3fi, while for SI interactions the XE10 bound still im- 
plies tension at more than 5tr for m x > 10 GeV. Hence, 
the precise CL of exclusion may depend on systematic 
uncertainties. 

In conclusion, we have presented a powerful method 
to check the consistency of an annual modulation signal 
in a DM direct detection experiment with bounds on the 
total DM scattering rate from other experiments, almost 
completely independent of astrophysics, for a given type 
of DM-nucleus interaction. While our bounds strongly 
disfavor a DM interpretation of present annually modu- 
lated signals in the case of SI and SD elastic scattering, 
the method will be an important test that any future 
modulated signal will have to pass before a DM interpre- 
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